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Abstract: We present realistic 3D numerical simulations of elastic bodies sliding 
on top of each other in a regime of velocities ranging from meters to tens of meters 
per second using the so-called Smoothed Particle Hydrodynamics (SPH) method. 
This allows us to probe intimately the response of the bodies and the nature of the 
friction between them. Our investigations are restricted to regimes of pressure and 
roughness where only elastic deformations occur between asperities at the contact 
surface between the slider block and the substrate. In this regime, solid friction is 
due to the generation of vibrational radiations which are subsequently escaping to 
infinity or damped out in which case energy is dissipated. We study periodic com- 
mensurate and incommensurate asperities and various types of disordered surfaces. 
In the elastic regime studied here, we report the evidence of a transition from zero 
(or non-measurable /i < 0.001) friction to a finite friction as the normal pressure 
increases above about 10^ Pa. For larger normal pressures (up to 10^ -Po), we find a 
remarkably universal value for the friction coefficient /i ~ 0.06, which is independent 
of the internal dissipation strength over three order of magnitudes, and independent 
of the detailled nature of the slider block-substrate interactions. We find that dis- 
order may either decrease or increase ^ due to the competition between two effects: 
disorder detunes the coherent vibrations of the asperties that occur in the periodic 
case, leading to weaker acoustic radiation and thus weaker damping. On the other 
hand, large disorder leads to stronger vibration amplitudes at local asperities and thus 
stronger damping. Our simulations have confirmed the existence of jumps over steps 
or asperities of the slider blocks occurring at the largest velocities studied (10 m/s). 
These jumps lead to chaotic motions similar to the bouncing-ball problem. We find 
a velocity strengthening with a doubling of the friction coefficient as the velocity 
increases from 1 m/s to 10 m/s. This refiects the increasing strength of vibrational 
damping. 



1 Introduction 



Solid friction has a long scientific history, starting probably in the western world 
about 500 years ago with the work of Leonardo de Vinci, continuing with the empir- 
ical Amontons'laws two centuries later and Coulomb's investigations of the influence 
of slipping velocity on friction in the XVIII century. Only three decades ago was it 
recognized that friction plays probably a fundamental role in the mechanics of earth- 
quakes [Brace and Byerlee, 1966]. Rock mechanicians consider an earthquake as a 
stick-slip event controlled by the friction properties of the fault, i.e. the destabiliza- 
tion of a weak part of the crust. This formulation has been shaped by laboratory 
experiments performed under a variety of pressure and temperature conditions (which 
however reproduce only imperfectly the conditions prevailing in the crust). Numerous 
laboratory experiments have been carried out to identify the parameters that control 
sohd friction and its stick-shp behavior [Persson and Tosatti, 1996]. The most signif- 
icant variables appear to be the mineralogy, the porosity, the thickness of the gouge, 
the effective pressure, the temperature and the water content [Byerlee et al, 1968; 
Brace, 1972; Beeman et al., 1988; Gu and Wong, 1991; Johansen et al, 1993; Streit, 
1997]. Low velocity experiments have established that solid friction is a function of 
both the velocity of sliding and of one or several state parameters, roughly quantify- 
ing the true surface of contact [Brace, 1972; Dieterich, 1972; 1978; 1979; 1992; Ruina, 
1983; Cox, 1990; Beeler et al, 1994; 1996; Baumberger and Gauthier, 1996; Scholz, 
1998]. 

The Ruina-Dieterich laws constitute the basic ingredients in most models and 
numerical elastodynamic calculations that attempt to understand the characteristics 
of earthquake sources. A recent lively debate has been whether space-time complexity 
in earthquake sequences can occur on an homogeneous fault solely from the nonlinear 
dynamics [Shaw, 1993; 1995; 1997; Gochard and Madariaga, 1994; 1996] associated 
with the slip and the velocity dependent friction law [Dieterich and Kilgore, 1994; 
Dieterich, 1992], or, does it necessarily require the presence of quenched heterogeneity 
[Rice, 1993; Benzion and Rice, 1993; 1995; Knopojf, 1996]? It is now understood that 
complexity can emerge purely from the nonlinear laws but heterogeneity is probably 
the most important factor dominating the multi-scale complex nature of earthquakes 
and faulting [Ouillon et al, 1996]. It is also known to control the appearance of 
self-organized critical behavior in a class of models relevant to the crust [Sornette et 
al, 1995; Shnirman and Blanter, 1998]. 

A well-known and serious limitation of these calculations based on laboratory fric- 
tion experiments is that the friction laws have been determined using sliding velocities 
no more than about 1 cm/s, i.e. orders of magnitude below the sliding velocity of 
meters or tens of meters believed to occur during earthquakes. The validity of extrap- 
olations, especially the velocity weakening dependence, has yet to be demonstrated. 
This is all the more relevant when one examines the underlying physical mechanisms 
giving rise to the friction laws. At low velocity, hysteretic elastic and plastic defor- 
mations at the length scale corresponding to asperities seem to play a dominant role 
[Bowden and Tabor, 1954; Sokoloff, 1984; Jensen et al., 1993; Dieterich and Kilgore, 
1994; Caroli and Nozieres, 1996; Tanguy and Nozieres, 1996; Tanguy and Roux, 1997; 
Garoli and Velicky, 1997; Bocquet and Jensen, 1997]. 

At larger velocities, different mechanisms come into play. Collisions between as- 



perities and transfer of momentum between the directions parallel and perpendicular 

to the motion may become an important mechanism [Lomnitz, 1991; Pisarenko and 
Mora, 1994]. This regime is probably relevant to explain the apparent low heat flow 
and thus low friction coefficient observed along the San Andreas fault, the so-called 
heat flow paradox. This paradox [BuUard, 1954] comes about because, in order to al- 
low for large earthquakes, a fault should have a large friction coefficient so that it can 
store large amount of elastic energy. However, repeated earthquakes occurring with a 
large coefficient of dynamical friction should give rise to a signiflcant heat flow at the 
surface which has not been observed [Henyey and Wasserburg, 1971; Lachenbruch and 
Sass, 1980]. One explanation for this low heat flow is that the coefficient of friction 
is low as a result of dynamical effects affecting the friction during the earthquake. 
Other possibilities, less explored in the literature, involve fluids [National Research 
Council, 1990] or acoustic vibrations [Melosh, 1996]. Several simplifled models have 
recently been explored as possible mechanisms for the generation of a low friction. 
These mechanisms include crack opening modes of slip [Brune et al., 1993], dynamical 
collision effects [Lomnitz, 1991; Pisarenko and Mora, 1994], frictional properties of a 
granular gouge model under large slip [Scott, 1996], space fllling bearings with com- 
patible kinematic rotations [Herrmann et al., 1990], hierarchical scaling [Schmittbuhl 
et al, 1996]. 

Recently, Tsutsumi and Shimamoto [1996; 1997] have reached a completely novel 
regime, by performing friction measurements on rotation cylindrical samples at ve- 
locities up to 1.8 m/s and for slips of several tens of meters. While these results due 
to several experimental problems are not completely straight-forward to interpret, 
they seem to indicate the existence of a change of regime from velocity weakening to 
velocity strengthening and then again to velocity weakening at the largest velocities. 
This last regime seems to be associated to the melting of a very thin layer. 

The work we report here has been motivated by considerations that different 
physical mechanisms might lead to a change of regime in the velocity and slip depen- 
dence of the solid friction law. Thus, one needs to explore the high velocity regime 
in as large a variety of conditions as has been done for the low velocity regime. For 
this purpose, we have developed realistic 3D numerical simulations of elastic bodies 
sliding on top of each other in a regime of velocities ranging from meters to tens of 
meters per second. In this way, we probe more intimately than any exprimental setup 
could do the response of the bodies and the nature of the friction. 

We begin by a short presentation of the SPH method and its implementation in 
our context. We follow by a presentation and a discussion of the results obtained 
for surfaces with periodic asperities and then with various types of random asperities 
and conclude. 

2 The numerical model 
2.1 The SPH method 

We have adapted the Lagrangian method called "smoothed particle hydrodynamics" 
(SPH), initially introduced by Lucy [1977] for hydrodynamic problems with fast dy- 
namics. This method has the avantage of being simple, elegant, easy to implement 



and to extend with a reasonable precision. Its most recent improved version makes 
it well-suited to treat problems with fast moving interfaces. This was our initial mo- 
tivation to use this property at our advantage to tackle the solid friction problem in 
the regime of large slipping velocities. 

In SPH the physical problem is discretized on a mesh whose nodes, the "particles" 
are moving or adapting in a Lagrangian manner. Each particles carries with it a set 
of field variables. An interpolating kernel allows to reconstruct the field variables 
everywhere by interpolation. Spatial derivatives are obtained from the analytical 
differentiation of the interpolation kernel. SPH has been used in a variety of appli- 
cations, such as gas dynamics [Monaghan and Gingold, 1983], fragmentation of gas 
clouds [Lattanzio et ai, 1985], radio jets [Coleman and Bicknell, 1985], impacts [Benz 
et ai, 1986; Benz and Asphaug, 1994; Benz et ai, 1994; Asphaug et ai, 1998], quasi- 
incompressible fluid flows [Monaghan and Humble, 1993], material rupture [Benz and 
Asphaug, 1995], dam rupture, ocean wave propagation and water falls [Monaghan, 
1994; 1996]. For reviews, see [Monaghan, 1988; 1992; Benz, 1990]. 

2.2 Outline of the approach 

Consider one of the physical field f{r), which is a function of position r. It can always 
be written as 

/(r1 = J f{r')S{r-r')dr' , (1) 

where S is the Dirac function. This suggests approximating f{r) by a smoothing 
kernel as 

= / fin W{r- r', h) dr' , (2) 

W{r — f',h) is the interpolating kernel and h is the width of the kernel and thus 
the smoothing scale. The kernel is continuous, differentiable and has the following 
properties: 

jw{f-f',h)dr' = l (3) 

and 

\im W {f - f' ,h) = 6{f - f') . (4) 

h — ^0 

/^From (^) and (^, (/(r)) ^^-^ /(^- From now on, we drop the arrow and write r as 
r, and similarly for the derivatives. For a fluid of density p(r), equation (j^) reads 

W{r -r',h)p{r')dr' . (5) 

Discretizing space in elements of masses mj, the integral becomes 

N e 

{f{r,)) = Y.m/-^W{n-r,,h) , (6) 
j=i Pj 
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with fj = f{rj). Replacing /(r) by p(r), this yields the following expression of the 
fluid density 



TV 



{p{ri)) = Y.m,W{r,-r,,h) . (7) 
i=i 

This equation has the following interpretation [Benz^ 1990] : each particle of mass 
ruj is spatially smoothed out according to the kernel space dependence, which can be 
seen as its spatial density distribution. The density at any point in space is obtained 
by summing all contributions from all particles at this point. The term Smoothed 
Particle Hydrodynamics expresses this interpretation. 
The gradient of / is similarly obtained as 



(V/(r)) = j Vf{r') W{r - r\ h) dr' . 



Integrating by parts and using the fact that W goes to zero sufficiently fast so that 
the surface terms are negligible, one obtains 



(V/(r)) = J VW{r - r', h) f{r') dr' , (9) 
with VW{r — r', h) being the gradient with respect to r. Discretizing, this yields 

(V/(r.)) = E— , (10) 
j=i Pj 

where Vj is the gradient with respect to the coordinate of the i-th particle and 
W,, = W{r,-r„h). 

Several choices are possible for the kernel, as long as the conditions (0) and (^) are 
fulfilled. Kernels constructed on spline functions have several advantages [Monaghan 
and Lattanzio, 1985] : 

^ I 1 - + if < < 1 
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else. 



where v = j^. This kernel has compact support, i.e. interactions between particles 
vanish for r > 2h. Only contributions from neighbors need to be accounted for 
instead of A^^ contributions. The second derivative of W{r, h) is continuous and the 
error in the estimation of the interpolation is of order 0{h'^). At the beginning of 
the calculation, the values of the kernel and its gradient for different values of v are 
calculated and stored in a array. This allows to decrease the calculation time. 



2.3 Equations of motion 

The formulation we use is borrowed from [Stellingwerf and Wingate, 1992; Benz and 
Asphaug, 1995]. The equation of mass conservation 



where v is the local velocity along a, gives after SPH discretization 

^-j:rn,{vr-vf)VtW,,. (13) 
The equation of momemtum conservation (neglecting gravity) 

where cr"'^ is the stress tensor defined by 

where P is the pressure, S"'^ is the deviatoric stress tensor with zero trace, and 5"^^ 
is the Kronecker symbol, becomes 

^=g™4^+^)vf«'«. (16) 

written in a symmetric form. It is easy to check that the total momemtum is conserved 
by multiplying this equation by rrii and verifying the exact symmetry in i and j. 
The equation of energy conservation 

dt p dx'^ p ^ ' 

where is the tensor of the rate of deformations defined by 

1 (dv'^ dvl^^ 



2 V dx^ dx'^ 

becomes 



+ (18) 



^^E"..K-<)(^)vfw'«. (19) 

This symmetric form ensures an exact energy conservation. 
The discretized expression of the deformation rate tensor is 

1 N 

^f= ^ E ^ (« - <) vf + « - v^mw^;) . (20) 



Thus 



j Pj dxi 

i (sr^rrij dWij rrij dWij , 

= 2 S -("J - + E -("I -r)^ I . (21) 



2.4 Deformation model 



The model of mechanical deformation that we implement is the simplest possible, 
namely a perfect elastic-plastic model, obeying Hooke's law in the elastic regime and 
a suitable plasticity criterion. In the results presented in this paper, we will not reach 
the regime where plasticity appears and leave the study of this regime for future work. 
In the elastic regime, the stress tensor reads 

^ = 2/i(e°'' --^"'^ e^^) + + S^^R''^ (22) 

at 3 

where /i is the shear modulus of the material and R is the tensor of stress rotations 
defined by 

^ ~ 2[dx(^ dx-) ■ ^^"^^ 
Its discrete SPH approximation is similar to that of the deformation rate tensor 

1 ^ m 

Rf = o E ((^; - <) Vf - {vl - v^mW,,) . (24) 

We need in addition to specify the equation of state, namely the dependence of 
the pressure P — P{p, u) as a function of the density p and the internal energy u. 
We use the equation of state due to Tillotson [Tillotson, 1962; Benz et ai, 1994], 
which works both for expanded as well as condensed phases under large pressure or 
impacts. This brings in the possibility to study the effect of melting or other extreme 
conditions that could occur locally under conditions of fast slipping rates. We will 
not fully exploit this possibility in the preliminary results presented below. We do 
find a local temperature rise at the level of boundary particles under friction but 
which is insufficient to lead to melting. We thus do not incorporate the physics of 
thermal diffusion and the effect of the internal energy is negligible in the regime of 
our simulations. The parameters used in our simulations are obtained from [Asphaug 
and Melosh, 1993] for typical geological rocks, such as granite, basalt and sandstone. 
Their densities are respectively 2.7 g cm^^ for granite and basalt and 2.3 g cm~^ for 
sandstone. We have not observed significant differences in the solid friction for these 
different materials. 



2.5 Geometry and scaling of the numerical experiments 

We consider the classical friction experiment in which a block of mass M in contact 
with a solid substrate is submitted to a normal pressure P and to a constant horizontal 
velocity v (see figure 1). We have worked with a block of size 0.5 x 0.375 x 0.25 cm^ 
while the substrate has dimension 1 x 0.5 x 0.25 cm^. Three codes are dedicated to the 
construction and positionning of the block and substrate. The first code constructs 
the objects as ensembles of particles. The second code chisels the boundaries of the 
block and substrate, their rugosities and put them in contact. The initial distance 
and conditions are thus determined. The third code calculates the final dimensions of 
the objects and retrieves the coordinates of the particles at the boundaries (first and 



The particles are put on a regular lattice. We have used the cubic lattice and the 
compact hexagonal lattice. In our simulations, we have used between 3000 and 15000 
particles, as a compromise between meshing and computation time. Figure 2 shows a 
configuration where the particles making up the block and substrate are represented : 
in this example, the total number of particles in the block is 12 x 16 x 20 = 1520. Thus 
the size of a particle is of order 0.025 cm. A slider block of a centimeter scale can thus 
be viewed as being constituted of effective particles of a fraction of a millimeter. We 
must thus incorporate the correct physics at the scale of each particle. Each particle 
at the boundary can be viewed as an elementary asperity that will interact with the 
particle-asperities of the substrate. 

Associated to their size h , mass m and interaction with neighbors, each particle 
has a characteristic oscillation frequency ujq given by 

^o = \ - ^T J- , 25 
V m V P 

where E is the Young modulus and p is the density of the material. We have used the 
fact that the effective elastic constant K felt by a particle is K ^ h E and the mass 
of a particle is m ^ p h^. Take h ~ 0.025 cm, E = 10^^ Pa and p ~ 3 10^ kg.cm~^, 
we get ti^o ~ 2.3 10^ and a natural frequency ujo/27t ^ 4 10® s^^. The expression 
(|25| ) can be written in terms of the period of oscillation 

^ 2-71 , nr h ^ , , 

To = — = 271 hJ— = 271 2.5 10"^ s , 26 
UJq \ E c 

where c is the longitudinal sound velocity. This period can be compared to the time 

h , , 

to = - 27 

V 

it takes for a particle driven at a velocity v to move over the distance h. Iq/Tq = 
(l/27r) {v/c) is thus small for subsonic slider block velocities. 

The boundary particles are accelerated due to their collision with the substrate 
under the imposed sliding velocity. Due to their acceleration, they entrain their neigh- 
bors which themselves accelerate and entrain their neighbors, and so on. Macroscop- 
ically, this is nothing but wave radiation which may act as an important damping 
process. Since the only physical ingredient of our model incorporates elastic interac- 
tion, friction can only emerge as a result of dissipation due to wave radiation. We thus 
need to estimate how the efficiency of this radiation is modified by the coarse-graining 
at the scale of the particles and correct for it. 

Recall that a particle of mass m submitted to a force Ff.xt accelerates according 

to 

rrii v= Fext ■ (28) 

Radiation is a generic phenomenon that refiects the acceleration of the body un- 
der consideration. Generically, the power radiated from a body accelerating at v is 
proportional to the square of its acceleration 



Pit) = mTn 



(29) 



Note that this resuh holds not only for acoustic waves, but for any wave (electro- 
magnetic, gravitational, hydrodynamics, etc). This radiation exerts a feedback force 
Frad that modifies the acceleration of the body as follows [Jackson, 1962]. Assuming 
the existence of this radiation force, we replace (pHj) by 



m v= Fe^t + Frad ■ (30) 

Frad is determined from the condition that its work during ti < t < t2 is equal to the 
radiation energy 

' Frad ■vdt = - r mTo dt . (31) 

tl Jti 



Integrating by part, we get 

rt2 rt2 



/ Frad ■ vdt = mTo / V ■ V dt — mTo{v -v) 



t2 

(32) 



For a periodic motion or if (?) -v) = at t = ti and t = t2, 

r {Frad - mTo v) ■ vdt = , (33) 
Jti 

thus leading to 

Frad = mTo d . (34) 

This effective force is indeed a dissipation as its sign changes under time reversal 
t ^ —t: recall that dissipation is nothing but the lack of invariance of the motion 
with respect to the change t —>■ —t. The important point is that the dissipation force 
due to radiation is proportional to the derivative of the acceleration, i.e. to the third- 
order derivative of the position. This is quite different from the first-order derivative 
dependence of standard fluid friction. The upshot is that radiation is extremely 
efficient at high frequencies, since its power is proportional to v and thus to the fourth 
power of the frequency, according to (^) (this is nothing but Rayleigh's scattering law 
for radiations from objects smaller than the wavelength, which is indeed universal as it 
relies solely on dimensional considerations as first derived by Rayleigh (see [Sornette, 
1989] for a review and references therein). The problem however from our perspective 
is that radiation efficiency becomes very small for small frequencies. Since coarse- 
graining using a finite particle size decreases the natural oscillation frequencies, the 
resulting radiated power will be largely reduced compared to the case of a real material 
in which the particle sizes are the atomic scale. 

To estimate this effect of coarse-graining on the radiation efficiency, we take v= 
AuJq with the typical amplitude of motion given hy A = h e. A reasonable estimation 
of the strain e is such that the elastic potential energy stored per unit volume (l/2)cre 



be equal to the kinetic energy density {l/2)pv'^. This leads to e = f \J p/E and thus 

A = h V \J p/E. Inserting in (^), we find that the radiated power per particle is 
P = 27T ^ p E h"^. Physically, the important quantity to establish the energy 
balance and the friction law is the radiated energy per unit volume 

2 



Since h is about 10^ larger than the atomic scale, the radiation efficiency in our 
coarse-grained formulation stemming from the acceleration of the SPH particles may 
be down to 10~^ that of expected physical one. This implies that we need to add a 
dissipation term to the equation in each particle to account for the physical radiation 
due to the accelerations at the sub-particle scale that is not correctly described by the 
coarse-graining at the particle scale. Our hypothesis, that will be conffimed below, 
is that friction is not sensitive to the specific form of the internal dissipation as long 
as it is present to damp out the vibrations. 

2.6 Dissipation 

Friction is fundamentally about dissipation. As we already mentionned, this dissipa- 
tion can be in the form of radiated waves or conversion of local vibrations (phonons) 
into others modes, thus corresponding to an effective loss of energy. To be consistent 
with our preceding discussion, we need to add a dissipation term in the equation of 
motion proportional to the derivative of the acceleration. We would however need to 
scale the amplitude of this term by a factor up to 10^ in order to obtain a correct 
scaling of the radiation power at the characteristic particle frequency. The problem 
however is that a much larger spectrum of frequencies are excited in the complex 
sliding motion and it is not possible to scale simultaneously the radiation at all these 
frequencies simultaneously. We have prefered a simpler approach which is to add a 
standard viscous dissipation force on each particle 

fdiss = -m'jv , (36) 

with a viscous coefficient which is a parameter of the model. Here, v is the velocity of 
a particle with respect to the center of mass of the slider block. As a consequence, this 
viscous dissipation acting within each particle damps their motion and thus converts a 
part of the sliding kinetic energy into losses that will finally produce the solid friction 
behavior. We have varied 7 in the interval between and 10^° s~^. Consistent with 
the use of this term (|36D as a device to mimick losses from radiation processes, we 
have not included this term in the energy balance equation. In other words, any heat 
produced by this dissipation is assumed to be instantaneously radiated. 

Before resorting to this viscous dissipation, we have explored various other possi- 
bilities, such as trapping the acoustic waves in dissipative cavities of various form, so 
as to mimick out-going escaping radiated waves that never come back. Unfortunately, 
in addition to the partial backscattering occuring at junctions (even with our best ef- 
fort to adapt the acoustic impedance by using smooth geometries and slowly varying 
mechanical properties), we found that the pattern of the acoustic particle vibrations 
self-organized so that the amplitude became vanishing small at the borders of these 
traps, making them totally inefficient. Recall that during a typical experiment of one 
millisecond at 10 m/s, the slider block slips over 1 cm while the acoustic vibrations 
propagate over about 5 m, i.e. have time to make 5000 travels back and forth within 
the slider block. There is a lot of shaking and organization of particle vibration going 
on all the time! 

From a numerical point of view, it is also necessary to incorporate a numerical 
viscosity that allows to regularize large gradients as occur for instance in shocks. This 



artificial viscosity modifies the equation of momentum conservation into 



^ = E-.Rr + ^ + nJvfw^,.. (37) 

j=l \ Hi Pj J 

Correspondingly, the artificial viscosity also brings in a contribution to the equation 
of energy conservation. The addition term is similar to a pressure 

n, = |"^==^ (38) 

[ else , 

where Cij = |(cj + Cj) and c, and Cj are the sound velocities in the particles i and j. 
The average density is defined by p^^- = |(pi + pj)- The estimation of the divergence 
of the relative velocity between i et j is 

h(vj — Vi) ■ ifj — fi) 

The addition term eh"^ prevents a divergence occurring for small \ fi — rj\. For instance, 
the choice e = 0.01 gives a correct smoothing of the velocity only if the interparticle 
distance is larger than 0.1/i. 

Note that liij remains symmetric in i and j which ensures the conservation of 
the total linear and angular momenta. The artificial viscosity disappears when the 
two particles separate, which ensures that this dissipation obeys the second law of 
thermodynamics, i.e. can only increase the entropy of the system. The values of 
the parameters a and (3 are not critical. In the following, we adopt a = 1 and 
(3 = 0.1. We have checked that the additional "dissipation" stemming from rounding 
and integrating errors is negligible. 

This artificial viscosity introduced in the SPH formulation may produce a non- 
negligible dissipation when shear is important and we have been concerned with the 
possibility that the effects reported below could stem purely from numerical effects. 
We believe that the artificial viscosity does not produce in itself a solid friction from 
three observations : 

• the relative velocities of the particles with respect to the center of mass are 
small ; 

• the solid friction coefficient is zero, to within numerical accuracy, when 7 = 0, 
i.e. the artificial viscosity alone does not produce a detectable friction ; 

• the temperature does not increase appreciably. This provides an upperbound 
for the dissipated energy due to the artificial viscosity (recall that the viscosity 
(^) is not incorporated in the heat production as mentionned above), which is 
thus found negligible. 



We conclude that the results reported below are not the expression of numerical 
tuning. 



2.7 Implementation 



Each particle carries 13 variables : three positions {x, y, z), three velocities {vx, Vy, Vz), 
the density, the energy, and five components of the deviatoric stress tensor Sxx, Syy, Sxy, Sxz, s. 
At each time step, this set of 13 variables is incremented. The temporal integration 
is carried out using a second order Runge-Kunta-Fehlberg algorithm with adaptative 
time steps. To ensure second order accuracy, the forces arc evaluated twice per time 
step. When a particle is suddenly colliding with another one, the local stress can 
jump to high values resulting in a drastic reduction of the time step reaching in some 
cases a factor 100! The typical time step is 5 10~^ seconds. Our simulations are run 
over 20, 000 to 50, 000 time steps, i.e. over a total duration of a millisecond or more. 
The simulations have been performed on RISK6000 and Ultra-Sparc 2 workstations 
and typically last three days. 

Starting from a configuration where the block is in contact with the substrate, 
at a distance equal to the inter-particle spacing, we submit the block to a vertical 
normal pressure along the z axis. This is done by applying the same force on all the 
particles on the upper boundary of the block. The bottom boundary of the substrate 
is fixed and cannot move. The pressure is applied progressively as 1 — exp(— t/r), with 
a characteristic r fsi 30 — 250 fxs, which is not too short to allow for the propagation 
of acoustic waves back and for the system to reach mechanical equilibrium. We have 
varied the pressure from 10^ Pa to 10^ Pa. Once the equilibrium is — bf reached, we 
impose a horizontal sliding velocity v. This velocity was varied v in the range 0.1 to 
10 m/s. During the first integration time step, all particle in the sliding block are 
moved with velocity v. During the rest of the integration, only those particles in the 
first and sometimes second layers of the sliding block are moved with velocity v the 
others ones evolving according to the laws of elasticity. This somewhat cumbersome 
starting procedure prevents large unwanted inertial oscillations to occur that would 
require a long time to damp and brings the slider block efficiently to a constant 
velocity. Our choice to push the sliding block at its upper boundary is an attempt 
to mimick large scale driving boundary conditions. It also allows the block to freely 
adjust itself close to the interface with the substrate in the hope to minimize the 
influence of boundary conditions. 

At each time step, the forces exerted on each particle are calculated. They are 
the cohesion force f^oh between each particle, the uniaxial pressure fpi, exerted at the 
upper boundary of the slider block, the repulsion force /^p,. at the block-substrate 
interface that acts only between the particles in the first layer of each objects, the 
viscous dissipation force fdis — —j{vx — v) along the x-axis, and —jVy along the 
y-axis and —jVz along the 2;-axis. The conservation of momentum reads 



The solid friction force is measured as the force /fp,,. exerted along the direction of 
motion Ox by the substrate on the particles in the first layer of the slider block in 
contact with the substrate. We also measure the vertical component and the ratio 



ITT' fcoh ~l~ fpb ~l~ fspr ~l~ fdis ■ 



(40) 



f spr / f ST 



spr 



(41) 



gives a local and instantaneous measure of the solid friction coefficient. We can then 



We have investigated several types of boundary forces between substrate and slider 
block particles : 



1. spring-like force 



Mr) = 



-K{r -h) <r < h 
r > h 



(42) 



2. Smooth repulsive force [Monaghan, 1994] 

/2(r) 



' A{^y Q<r<l 
D{r-hf l<r<h 
r>h 



(43) 



This force is also radial and its first derivative is continuous. The continuity 
condition r — I gives D — A{^Y and I — j^^- In our calculations, we 
take b — h. 



3. Lennard- Jones force per unit of mass 
Mr) = < 



e ((^)"' - (^)") < r < Tc with m>n 
A{R-rf + D{R-r) Tc < r < R 
r > R 



(44) 



Mr) vanishes for r = Z and r — R and has a minimum at r = Tc- In our 

simulations, we have taken R = |Ap and I — Ap. The continuity of Mr) ^nd 



of its first derivative at r = Tc implies 
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from which we derive 
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e remains an adjustable parameter that is chosen so that the time scale 

1 h 



With Ap = h , this leads to 
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(47) 
(48) 
(49) 



We have used m = 8,n — 4 and m — 12,n — 6 and found no significant 



Even without introducing asperities explicit, the potential field seen by particles 
at the interface between substrate and slider block is not smooth. Due to the particle 
structure, the potential thus exhibits a periodic structure of troughs and peaks with a 
period equal to that of the cubic or hexagonal lattice used to construct the block and 
substrate. This corresponds to a slider block with a periodically modulated roughness 
sliding on a substrate presenting also a periodically modulated roughness. We have 
also investigated cases for which we explicitly introduced random asperities. 



3 Solid friction for periodic roughness 
3.1 Measurements 

We have carried out simulations for various dissipation coefficient 7. Figure 3a shows 



the local friction coefficient n = fspr/fspr ^ defined in (^) measured on a single par- 



ticle in the first layer of the slider block in contact with the substrate. The simulation 



uses the boundary force /i defined by (42). The total number of particles used is 



5504, the sliding velocity v = 1 m/s, the applied pressure 10^ Pa, and the dissipation 
parameters a = 0.1, 7 = 10® s~^. The positive and negative oscillations correspond to 
alternative braking and accelerating phases of the particle as it climbs up and down 
the asperities of the substrate. The time average of this local instantaneous friction 
coefficient is /i ~ —0.05, the negative sign corresponding to a net drag. Figure 3b 
shows the global instantaneous friction fii, obtained by taking the ratio of the total 
force along x on all block particles on the boundary in contact with the substrate 
to the total force along z. We obtain the same estimate fii ~ —0.05 for the friction 
coefficient after time averaging. 

The following table summarizes our results for the effective friction coefficient 
for three different values of the dissipation 7. The simulations have otherwise been 
carried out exactly in the same way, using 5504 particles, an imposed sliding velocity 
of 1 m/s and an applied pressure of 10^ Pa. In the Table, 5/i denotes the standard 
deviation of /i, i.e. the amplitude of its fluctuations 



a 





0,1 


0,1 


7 








10® 


/i 


-0.020 


-0.016 


-0.058 


6fi 


0.017 


0.017 


0.007 




-0.021 


-0.031 


-0.050 



These results are compatible with a vanishing average friction in absence of internal 
dissipation. When internal friction is present, the erratic motions of the particles 
submitted to damped multiple acoustic paths decrease somewhat and produce a flnite 
friction coefficient which is significantly smaller than the instantaneous fiuctuations. 
Interestingly, the value of of the friction coefficient does not appear to change with 
increasing 7. Above 7 = 10^ s~^, the numerical time step becomes so small making 
the calculation almost impossible and thus preventing us from studying this regime. 
For 7 < 10® s~^, the fiuctuations are too great to allow for accurate measurements. 
To summarize, we obtain 

~ ,,. ~ _n nfi -1- n ni fr^^ 1 n6 „-i ^ ^ 1 n9 „-i ^'cn^ 



This result justifies a posteriori our procedure to model the radiation damping by 
the viscous dissipation (|36D . The remarkable fact that the solid friction coefficient is 
independent by and large of the amplitude of damping suggests the following picture: 
in the elastic regime we explore, solid friction stems from the acceleration of asper- 
ities brought in contact and collision that radiates high-frequency vibrational waves 
subsequently damped out, thus converting a part of the kinetic energy of the slider 
block into dissipation. The specific form of the internal dissipation seems not to be 
important, as long as there is a dissipation that can damp out the vibrations of the 
asperities. 



3.2 Particle motions 

Figure 4 shows the vertical motion of a particle of the slider block in the layer in 
contact with the substrate, using the boundary force /i, 5504 particles, v = 1 m/s, 
a pressure of 10^ Pa, a = 0.1,7 = 10^ s~^. The climb and fall over the particle- 
asperities of the substrate are clearly visible. The downward drift is caused by a 
steady horizontal drift of the slider block along the y direction, perpendicular to the 
driving velocity. The horizontal motion of such a typical particle is almost steadily 
increasing, with however a short stop just when the particle is at the bottom of the 
potential created by the substrate asperities. 

Figure 5 shows the velocity along x (fig. 5a) and along z (fig. 5b) of the center 
of mass of the slider block, under the same conditions as described in figure 4. This 
shows that the slider block as a whole moves up and down as well as periodically 
accelerates and decelerates due to the interactions with the periodic array of asperities 
of the substrate. 

For small 7 ~ 10® s~^, the slider block exhibits rather large vertical oscillations 
that decrease significantly in amplitude as 7 increases. For the largest explored 7, 
the slider block follows very closely the geometry of the substrate asperities. Note 
that the motion of a particle in the layer in contact with the substrate is a very good 
proxy for the motion of the slider block ClS db whole. 



3.3 Pressure and velocity dependence 

For pressure below 10'' Pa, we are unable to measure a non-zero friction coefficient. 
Independently of the dissipation level 7, the measured /i are within uncertainty the 
same with or without dissipation. The reason is that the asperities do not much 
penetrate into each other and the slider block "floats" over the substrate without 
generating significant vibrations that can be dissipated. While more simulations are 
required to demonstrate it, we surmise that this behavior is due to the existence in 
the elastic regime of a critical pressure threshold below which there is zero friction. 

For pressure above 10^ Pa, the substrate force /i cannot be used anymore as the 
slider block penetrates within the substrate. We have then used /2 and /s given by 



and (|4^) respectively. Unfortunately, the integration time step shrinks drasti- 
cally, thus limiting an exhaustive exploration of this regime. However, we have found 
that the results are the same for a pressure of 10^ Pa as found for 10^ Pa. We have 
also verified that the three different forces give the same results. 



With respect to the velocity dependence, our time-exphcit numerical SPH method 
docs not allow us to explore too small velocities due to the prohibitive calculation 
time. A few runs at w = 0.1 m/s give essentially the same results as for v = 1 m/s. 
However, for larger velocities v = 10 m/s, the friction coefficient increases and doubles 
at /i = 0.11. The instantaneous value fluctuates with much larger amplitudes as a 
result of a very high level of vibrational excitations induced by the collision between 
asperities. The corresponding instantaneous friction coefficient measured on all the 
particles in the boundary layer in contact with the substrate is shown in figure 6. The 
simulation uses the boundary force /i and has 5504 particles. The parameters are 
V — 10 m/s, which was imposed t — 0.2 ms, a pressure of 10^ Pa and dissipation 
q; = 0, 1, 7 = 10^ s~^. Note the existence of flat steps in the graph at the value ji — O: 
they correspond to jumps of the slider block over the substrate asperities. In these 
regimes, the slider block is literally flying over the substrate, as a result of an efficient 
transformation of horizontal to vertical momentum induced by the collisions with the 
substrate asperities. This regime has been postulated first by Lomnitz-Adler [1991] 
and our simulations confirm nicely his ideas. As a result of these jumps, the landing 
of the slider block does not occur in phase with the substrate. As a consequence, the 
evolution becomes chaotic, in the rigorous mathematical meaning of the term. The 
mechanism for this chaotic behavio is similar to that in the toy model of a bouncing 
ball on a sinusoidally vibrating table [Mehta and Luck, 1990; 1993; Franaszek and 
Isomaki, 1991; Luo and Han, 1996; de Oliveira and Goncalves, 1997]. The jump of 
the slider block occurs from roughly the maximum of the potential created by the 
substrate asperities and over its descent, i.e. in the pulling portion of the potential. 
This explains why the total time average friction coefficient is stronger as the pulling 
part as become weaker. We have not push more the numerical exploration of this 
very interesting behavior and leave it to a future publication. We expect even more 
interesting behavior at still larger velocities as the jump can carry the slider block 
over two or more asperities leading to the possibility of a rich phenomenology for the 
friction coefficient at these high sliding velocities. 

To summarize, the main result of our investigation of the velocity dependence of 
the solid friction coefficient is that it increases at large velocities. This is in agreement 
with the expectation that the vibrational radiation damping becomes the dominating 
mechanism with an efficiency that increases fast with the velocity. 

4 Disordered and fractal interfaces 

We have investigated three types of disordered roughness : incommensurate periodic 
roughness between the block and substrate, a step and random roughnesses. 

4.1 Commensurate periodic roughness 

Nothing changes compared to the previous periodic case, except for the fact that we 
tilt the lattice structure of the slider block with respect to the substrate by an angle 
between (previous case) and 45 degrees. When the angle is non zero, the asperities 
of the slider block do not encounter those of the substrate in the same configuration 
and at the same time. We find that, for most angles, the friction coefficient is slightly 



decreased compared to the value 0.06 of the periodic case. We find a remarkable 
result only for the special case of the most incommensurate regime where the ratio of 
the mesh size of the slider block to that of the substrate lattice is equal to the golden 
mean (\/5+ l)/2 ^ 1.618. Recall that the golden mean is the irrational number that 
is the least well approximated by a rational number. For this ratio, no two asperities 
of the slider block will be in the same configuration at the same time with respect 
to an asperity of the substrate. The measured friction coefficient ai v — 1 m/s, 
a normal pressure of 10® Pa and 7 = 10® s^^ is extremely small : /// = —0.003. 
Qualitatively, we attribute this small value to the conjunction of two effects. First, 
the tuning of the vibrational resonances occurring in the perfect periodic case does 
not occur anymore. This leads to much smaller coherent vibrations and thus smaller 
damping. Secondly, the still regular smooth roughness does not produce large local 
amplitudes of vibrations. We are nevertheless surprised that these effects contribute 
to such a small value of /ii. Our result is reminiscent of the "super- lubrification" 
regime found recently [Hirano et al, 1997] using a tungsten ^^(oii) wire sliding over 
a silicium S'i(ooi) surface in an incommensurate geometry and also over solid M0S2 
sohd films [Martin et al, 1993]. 

4.2 Step-like roughness 

The substrate is made of a "smooth" periodic surface up to some fixed Xstep-, at which 
a vertical step equal to the particle size h is made by adding a single particle layer 
beyond Xgtep- The simulations are performed as before. The slider block is accelerated 
at V before the step. The simulations use the substrate force /i, 5328 particles forming 
cubic lattices. Figure 7 corresponds to v = 1 m/s and a normal pressure of 10^ Pa. 
The arrows represent the instantaneous velocities of the particle. Figure 7a shows 
a snapshot exactly at the time when the slider block encounters the step. Figure 
7b shows that the slider block is ejected vertically and starts to jump over the step. 
Figure 8 corresponds to v — 10 m/s and a normal pressure of 10^ Pa. It shows a 
latter stage when the slider block is in fiight above the substrate. The arrows now 
show the stress carried by each particle projected in the 2D {x, z) plane. As the slider 
block flies over the substrate, the stress within it relaxes to zero. Its landing occurs 
several particles after the step and the slider block is found to bounce back several 
time before resumming its steady state sliding. These simulations demonstrate again 
the importance of jumps at high velocities, even in the presence of strong confining 
pressure. 

4.3 Random roughnesses 

Three types of randomness have been investigated: holes in the first layer of the sub- 
trate, variable heights of the substrate particles in the first layer and fractal roughness. 

4.3.1 Holes in the first layer of the substrate 

One could imagine first to introduce disorder by removing at random a fraction of 
the particles in the first layer of the substrate in contact with the slider block, thus 
creating holes of varying sizes and shapes controlled by the distribution of cluster 



sizes in 2D percolation [Stauffer and Aharony, 1994]. It turns out that, for a density 
of holes no larger than 60 %, this has no effect on the slider block as it continues to be 
supported by the remaining particles of the substrate. We find the same coefficient 
of friction as in absence of holes. Above 60 %, the slider block starts falling partly 
in sufficiently large holes and the regime of sliding becomes controlled by the jumps 
over steps as just described. This value of 60 % corresponds approximately to the 
concentration of holes at which the large clusters in the substrate become of size 
comparable to the slider block. 

4.3.2 Variable heights of the substrate particles in the first layer 

The substrate is now made of a single layer of particles. These particles are again 
positionned regularly on a lattice in the x — y plane, but their vertical positions are 
taken randomly and uniformly between —Az/2 and +Az/2. The shder block is not 
modified. The solid friction coefficient is now measured by measuring the total force 
exerted on the first layer of the slider block in contact with the substrate, as the 
particles in the first layer of the slider block are not continuously in contact with 
the substrate due to its random roughness. Figure 9 shows the measured friction 
coefficient ni as a function of time for a simulation using /i, 6640 particles forming 
a cubic lattice, v = 1 m/s and a normal pressure of 10^ Pa. After an initial large 
resistance at the beginning of the slider block motion, the friction coefficient settles to 
a stationary regime characterized by random fluctuations still decorated by a periodic 
structure reflecting that of the asperities of the slider block. We flnd that ni increases 
with the roughness amplitude Az of the substrate. For Az — 0.5 h, jii — —0.017 and 
increases continuously to ni = —0.08 for Az — h. Surprisingly, a small roughness 
decreases the solid friction while a larger roughness increases it above its periodic 
roughness value 0.06. We attribute the decrease of fj,i for small Az to the detuning 
of the vibrational resonances occurring in the perfect periodic case, that were at the 
origin of relatively large vibrational radiation and thus damping. 

4.3.3 Fractal roughness 

As in the previous section, the substrate is made of a single layer of particles, po- 
sitionned regularly on a 2D cubic lattice in the x — y plane with mesh h. Their 
vertical positions are determined by using the spectral synthesis method described 
in [Peztgen and Saupe, 1988] to generate a self-affine surface. We have investigated 
different dimensions between Df — 2.1 to Df — 2.9. The shder block is made of 6440 
particles organized in a hexagonal compact lattice. Thus, even without the fractal 
vertical structure of the substrate, there is no more commensurability between the 
slider block and substrate. The maximum amplitude of the self-affine surface is im- 
posed equal to h for all values oi Df. The largest wavelength that we have kept in 
the construction of the substrate is equal to one eighth the length of the shder block. 
This ensures that the slider block remains stable and does not fall or jump over steps 
as in the step case. Figures 10a and 10b show two fractal surfaces, respectively with 
Df = 2.3 and Df = 2.8, with the same maximum amplitude h. Figures 11a and 
lib show the locii of contacts between the substrate and the shder block for the two 
surfaces shown in figure 10. Figure 11a, corresponding to Df — 2.3, shows a larger 



and more coherent contact area than figure lib, corresponding to Df — 2.8. 

Figure 12 shows the variation of the sohd friction coefficient function of 

Df for simulations performed under a normal pressure of 10*^ Pa, v = 1 m/s and 
7 = 10® s^^. For Df < 2.6, the friction coefficient is found less than that of the perfect 
periodic case. As for the previous case, we attribute this result to the detuning of 
the vibrational resonances occurring in the perfect periodic case, that were at the 
origin of relatively large vibrational radiation and thus damping. For Df > 2.6, the 
friction coefficient becomes larger than that of the perfect periodic case, reflecting 
the transition to another regime controlled by the larger elastic distorsion of the 
asperities. 

5 Conclusion 

We have presented a new numerical model for the investigation of solid friction prop- 
erties in the regime of fast relative velocities of the order of meters to tens of meters 
per second. We have restricted our investigation to the regime where only elastic 
deformations occur between the asperities at the contact between the shder block 
and the substrate. In this case, the only mechanism that dissipates energy and cre- 
ates a non- vanishing solid friction coefficient is through the generation of vibrational 
radiations that are subsequently damped out, either by escaping to infinity or by a 
suitable internal damping process. We have examined periodic commensurate and 
incommensurate asperities and various types of disordered surface. In this elastic 
regime, we report the evidence of a transition from zero (or non-measurable) friction 
to a finite friction when the normal pressure increases above about 10® Pa. We find a 
remarkably universal value for the friction coefficient fi ~ 0.06, which is independent 
of the internal dissipation strength over three order of magnitudes, and independent 
of the detained nature of the slider block-substrate interaction. We find that disor- 
der may either decrease or increase fi due to the competition between two effects: 
Disorder detunes the coherent vibrations of the asperties that occur in the periodic 
case, leading to weaker acoustic radiation and thus weaker damping. On the other 
hand, large disorder leads to stronger vibration amplitudes at local asperities and 
thus stronger damping. Our simulations have confirmed the existence of jumps of 
the slider blocks that occur over steps or asperities for the largest velocities of 10 m/s. 
We find a velocity strengthening with a doubling of the friction coefficient when the 
velocity increases from 1 m/s to 10 m/s. This reflects the increasing amplitude of 
vibrational damping. 

We leave for another a later paper the investigation of the regime where the 
pressure is larger and the roughness is more disordered so that the local pressure 
at asperities reaches the plasticity threshold. In this regime, temperature diffusion 
must be added to the formulation. This does not pose any conceptual difficulty and 
can easily be incorporated in our code. In this regime, both vibration damping and 
plasticity become the source of dissipation. It is probable that the friction coefficient 
will be found larger in this often measured in macroscopic experiments that 

work in the regime where asperities are deformed in the plastic regime [Dieterich and 
Kilgore, 1994]. However, we still expect that jumps play an important role at the 
largest velocities of tens of meter per second. 



We acknowledge stimulating discussions with T. Villemin and J.-P. Vila. 
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Figure captions 



Fig. 1: Classical friction experiment in which a block of mass M in contact with 
a solid substrate is submitted to a normal pressure P and to a constant horizontal 
velocity v. We have worked with a block of size 0.5 x 0.375 x 0.25 cm? while the 
substrate has dimension 1 x 0.5 x 0.25 cm?. 



Fig. 2: A configuration where the particles making up the block and substrate 
are represented. In this example, the total number of particles in the block is 12 x 
16 X 20 = 1520. The size of a particle is of order 0.025 cm. A slider block of a 
centimeter scale is constituted of effective particles of a fraction of a millimeter that 
act as constitutive grains. Each particle at the boundary is an elementary asperity 
that interacts with the particle-asperities of the substrate. 

Fig. 3: a) Local friction coefficient /i = fgpr/fspr defined in (^) measured on 
a single particle in the first layer of the slider block in contact with the substrate. 
The simulation uses the boundary force /i defined by (|42|). The total number of 
particles is 5504, the sliding velocity is f = 1 m/s, the applied pressure is 10® Pa, 
the dissipation parameters are a = 0.1,7 = 10^ s~^. The time average of this local 
instantaneous friction coefficient is yU ~ —0.05, the negative sign corresponding to a 
net drag. 

b) Global instantaneous friction /i;, obtained by taking the ratio of the total force 
along X on all block particles on the boundary in contact with the substrate to the 
total force along z. We obtain the same estimate fii ~ —0.05 for the friction coefficient 
when time averaging. 

Fig. 4: Vertical motion of a particle of the slider block in the layer in contact with 
the substrate, using the boundary force /i, 5504 particles, v = 1 m/s, a pressure of 
10® Pa, a = 0.1,7 = 10^ s-\ 

Fig. 5: Velocity along x (fig. 5a) and along z (fig. 5b) of the center of mass of the 
slider block, under the same conditions as for figure 4. 

Fig. 6: Instantaneous friction coefficient measured on all the particles in the 
boundary layer in contact with the substrate using the boundary force /i, 5504 parti- 
cles, f = 10 m/s which was imposed at t = 2 ms, a pressure of 10® Pa and dissipation 
a = 0, 1,7 = 10^ s~^. The fiat steps in the graph at the value = correspond to 
jumps of the slider block over the substrate asperities. 

Fig. 7: Motion of the slider block over a step of height h after having been 
accelerated v = 1 m/s before the step. The simulations use the substrate force /i, 
5328 particles forming cubic lattices and a normal pressure of 10^ Pa. The arrows 
represent the instantaneous velocities of the particle. 

a) Snapshot exactly at the time when the slider block encounters the step. 

b) The slider block is ejected vertically and starts to jump over the step. 
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V — 10 m/s and a normal pressure of 10^ Pa. The arrows show the stress carried by 
each particle projected in the 2D {x, z) plane. 

Fig. 9: Measured friction coefficient function of time for a simulation using 

/i, 6640 particles forming cubic lattices, v = 1 m/s and a normal pressure of 10^ Pa. 
The substrate is made of a single layer of particles positionned regularly on a lattice 
in the x — y plane, but with their vertical positions taken randomly and uniformly 
between —Az/2 and +Az/2. The slider block is not modified. 

Fig. 10: Two fractal surfaces, respectively with fractal dimensions Df — 2.3 
(fig. 10a) and Df — 2.8 (fig. 10b), with the same maximum amplitude equal to h. As 
for fig. 9, the substrate is made of a single layer of particles, positionned regularly 
on a 2D cubic lattice in the x — y plane with mesh h. Their vertical positions are 
determined by using the spectral synthesis method described in [Peitgen and Saupe, 
1988] to generate a self-affine surface. The slider block is made of 6440 particles 
organized in a hexagonal compact lattice. The largest wavelength in the fractal 
surfaces is equal to one eighth the length of the slider block to ensure the stability of 
the slider block during its motion. 

Fig. 11: Locii of contacts between the substrate and the slider block for the two 
surfaces shown in figure 10. Figure 11a, corresponding to Df = 2.3, shows a larger 
and more coherent contact area than figure lib, corresponding to Df — 2.8. 

Fig. 12: Variation of the sohd friction coefficient ni as a function of the substrate 
fraction dimension Df for simulations performed under a normal pressure of 10^ Pa, 

V — 1 m/s and 7 = 10^ s~^. 
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